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Abstract A formalism is presented that allows an asymptotically exact solution of 
non-relativistic and semi-relativistic two-body problems with infinitely rising confining 
potentials. We consider both linear and quadratic confinement. The additional short- 
range terms are expanded in a Coulomb-Sturmian basis. Such kinds of Hamiltonians 
are frequently used in atomic, nuclear, and particle physics. 
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1 Introduction 



Non-relativistic and semi-relativistic Hamiltonians with infinitely rising potentials are 
often used to model quark-quark interactions or other quantum particles in a trap. 
There are many methods for solving the Schrodinger equation with these kinds of 
potentials. They mostly rely on variational principles and can provide satisfactory 
results for the few lowest-lying states. However, they are unable to account for the 
most important property of the system; they cannot grasp the infinitely many bound 
states, which can be important if we want to consider three or more particles. 

In this paper we offer a method that can provide an exact treatment of the infinitely 
rising confinement potentials. We consider both non-relativistic and semi-relativistic 
kinetic energy with linear and quadratic confinement. We write the bound-state prob- 
lem in integral equation form and incorporate the kinetic energy and the confining 
terms into the Green's operator. Then we perform a separable expansion of the short- 
range terms in the Coulomb-Sturmian basis. This basis allows an exact evaluation of 
the Green's operator in terms of continued fractions and provides an asymptotically 
exact treatment of problems with linear and quadratic confinement. 
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This paper is organized as follows. In Sec.[2]we devise our method for non-relativistic 
kinetic energy with linear or quadratic confinement. Sec. [3] is devoted to the semi- 
relativistic case. In Sec. [4] we present some numerical illustrations. We consider both 
linear and quadratic confinement with non-relativistic and semi-relativistic kinetic en- 
ergy. Finally, in Sec. |5j we summarize our results. 



2 Non-relativistic kinetic energy with short-range plus confining potential 

We consider first a Hamiltonian 

h = h° + v (c) + v (s) , (I) 

where h° is the non-relativistic kinetic energy operator, and v^- 8 ' are confining and 
short-range potentials, respectively. The non-relativistic kinetic energy is given by 

h° = ^- . (2) 
2m v ' 

For we assume a polynomial potential with linear and a quadratic confinement 

JT C ' = a—i/r + ao + a\r + air 2 . (3) 

We will see that the 1/r and a constant term in does not make the evaluation of 
the Green's operator any harder. 

To set up an integral equation, we need the asymptotic Hamiltonian. In this case 
the asymptotic Hamiltonian can be identified as 

h c = h a + « (c) . (4) 

Then the original Hamiltonian becomes 

h = h c + v {s) . (5) 

The bound states are the solutions of the Lippmann-Schwinger equation 

\^)=g c (E)v^\^) , (6) 

where 

g c (z) = (z-h c r 1 ■ (7) 



2.1 The Coulomb-Sturmian separable expansion method 



We solve the integral equations by using the Coulomb-Sturmian (CS) separable expan- 
sion approach. In configuration space, the CS functions are defined by 

{r\nl; b) = JZ + ^ + ^ (2br) l+1 e X p(~br)Ll l+1 (2br) , (8) 

where n = 0, 1 . . ., I is the angular momentum, 6 is a parameter and L is the Laguerre 
polynomial. The biorthonormal partner to \nl;b) is defined as 

{r\rdfb) = l/r{r\nl;b) . (9) 
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The orthogonality and completeness of the basis states, in angular momentum sub- 
space, is given by 

(n'f; b\nJfb) = (n^T&K 6) = 5 n , n 6 ln (10) 

and 

N N 
lim V \nl;b)(nl;b\ = lim V \nl; b){nl; b\ = lim if = l z , (11) 

n=0 n=0 

respectively. The CS functions also have a nice analytic form in momentum space 
,,, M /2 n! (rz + / + l);!b(4bp)' +1 ;+1 / p 2 -b 2 \ 

and 

^ = 2b(n + lll) ^- (13) 

We can approximate the short-range potential in a separable way. Here we 
adopt the double-basis expansion technique proposed originally in Ref. [T] and extended 
later to Coulomb-like potentials in Ref. 2 . It amounts to expanding as 

u (s) « ^ \nl; b) (<nZ; b\ml; b')j (ml-b'\v {s) \m'l'-b')[{m'l'-b'\n'U-b)^ (n'l';b\ 

nmm'n' 
N 

« ^|n/;% W (n'/';6| . (14) 



This is an exact representation if JV goes to infinity and becomes an approximation 
if N is kept finite. The CS matrix elements of the potential have to be evaluated 
numerically. This can easily be done in configuration or momentum space, depending 
how the potential is defined. 

If we plug the approximated potential operator into the Lippmann-Schwinger equa- 
tion we get 

N 

\TP) = J29 C (E)\nl;b)v^{n'l';b\il,) . (15) 

nn' 

Acting with the bra (n"/"; b\ from the left, we arrive at a matrix equation for the vector 

± = 9 c (E)v^±, (16) 

where g c (E) = {nl; b\g c (E)\n'l; b). In fact, this equation is a homogeneous linear equa- 
tion 

(l-£ c (£k (s) )^ = , (17) 
whose solution can be found as zeros of the determinant 

\(g c (E))- 1 -v^\=0. (18) 
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2.2 Green's operator for Hamiltonians with confinement potential 

The Green's operator |7|) is defined by the relation 

J(z) g c (z) = l, 

where 

j, \ P 2 a-l 2 

J(z) = z — - air — a,2r . 

2m r 

This equation in the CS basis looks like 

J2(nl;b\J(z)\n'l;b) (n^b\g c {z)\nFT; b) = S nn „ 



(19) 
(20) 

(21) 



and we need the N x N upper left part of the oo x oo Green's matrix. 

The operator J has an infinite matrix representation. It can be constructed from 
the following matrix elements [4]: 



(nl; b\ — \ril; b) — {n'l; b\ — \nl; b) — 5 nn i 
(nl; b\n I; b) = (n'l; b\nl; b) 



l(n + l + l) for n' = n 

- i y/n'(n' + 21 + 1) for n' = n + 1 
1 for n' > n + 1 



(22) 
(23) 



2 2 

(nl;b\- — \nl;b) — (nl;b\- — \nl;b) 
2m 2m 



'£r l ( n + l + 1 ) for n' = n 

-^^Jri(n' + 21 + 1) forn' = n+l (24) 
I for ri > n + 1 



{n'l; b\r\nl; b) = (nl;b\r\n l;b) 

b 
2 



b (n + l + l) 



for n' — n 
for n' — n + 1 



-\^n'{n' + 21 + 1) 
jij y/n'(n' - l)(n> + 21) (n' + 21 + 1) for n' = n + 2 
for ri > n + 2 



(25) 



{n'l; b\r 2 \nl; b) = (nl; b\r 2 \n'l; b) = 

' - [(((10n + 21 + 4) (n + 21 + 3) + 

9n(n - l))(n + 21 + 2) + n(n - l)(n - 2))] 



= < 



fff[(4n' + 2l)(ri + 21 + 2) + (ri - l)(n - 2)}^/n'(n' + 21 + 1) for n 
Vn'(n' - l)(n' + 21) (n' + 21 + 1) 
Is (2ri + 21) sjn'(n' - l)(n' - 2)(n' + 21 + l)(n' + 21) 



for 


n 


= n 


for 


n 


= 71+1 


for 


n 


= n + 2 


for 


n 


= n + 3 


for 


n 


> 71 + 3 






(26) 



In our previous study we found that the inverse of an infinite symmetric tridiagonal 
matrix can be determined by a continued fraction [3]. In this case, if the r 2 term 
is present in J, its matrix representation is a septadiagonal infinite symmetric band 
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matrix. Such a matrix can be considered as a tridiagonal matrix of 3 x 3 block matrices. 
If r is the highest power in J, the matrix is pentadiagonal, which can be considered as 
tridiagonal matrix of 2 x 2 blocks. An infinite symmetric band matrix can be considered 
as a tridiagonal matrix of block matrices. The corresponding Green's matrix can be 
constructed by matrix continued fractions [4]: 

!l( z ) = (J — &i,N'&i,N'JN' ,N' + lCN' + lJN' + l,N') j (27) 
where C is a matrix continued fraction defined recursively by 

Cj+1 = {Ji+l,i+l — Ji+l,i+2Ci+2Ji+2,i+l) 1 (28) 

and Ji j refers to 3 x 3 or 2 x 2 blocks. The index TV' = 7V/3 or TV' = N/2 depending 
on whether we cover J by 3x3 or 2x2 blocks. We found, that g c is an inverse of 
the modified J matrix. The modification is a matrix continued fraction C, which only 
effects the bottom right 3 x 3 or 2 x 2 block of J. 



3 Semi-relativistic kinetic energy with short-range plus confining potential 

Now we consider the problem of the previous section with semi-relativistic kinetic 
energy 

h = h° r + t/ c) + w (s) , (29) 

where 

h" = x/m 2 + p 2 - m . (30) 

The semi-relativistic kinetic energy operator in the CS basis does not have a tridiagonal 
or band-matrix structure. Numerical studies show that its matrix representation for 
large n and nf indices is pentadiagonal dominant (see Table [l]). These asymptotically 
dominant elements can easily be taken into account in our matrix continued fraction 
method formalism. 

Table 1 The CS matrix elements of h® = \/m 2 + p 2 — m for n,n' = 170 — 179. The penta- 
diaognal elements are clearly dominant. We used m = 0, I = 0, and b = 4. 



80. 


10. 


-30. 


-4. 


-5. 


-2. 


-2. 


-1. 


-1. 


-1. 


10. 


81. 


10. 


-30. 


-4. 


-5. 


-2. 


-2. 


-1. 


-1. 


-30. 


10. 


81. 


10. 


-30. 


-4. 


-5. 


-2. 


-2. 


-1. 


-4. 


-30. 


10. 


82. 


10. 


-30. 


-4. 


-5. 


-2. 


-2. 


-5. 


-4. 


-30. 


10. 


82. 


10. 


-31. 


-4. 


-5. 


-2. 


-2. 


-5. 


-4. 


-30. 


10. 


82. 


10. 


-31. 


-1. 


-5. 


-2. 


-2. 


-5. 


-4. 


-31. 


10. 


83. 


10. 


-31. 


-4. 


-1. 


-2. 


-2. 


-5. 


-4. 


-31. 


10. 


83. 


10. 


-31. 


-1. 


-1. 


-2. 


-2. 


-5. 


-4. 


-31. 


10. 


84. 


10. 


-1. 


-1. 


-1. 


-2. 


-2. 


-5. 


-4. 


-31. 


10. 


84. 



6 



4- 



2 - 







-2 - 



-4- 




Fig. 1 Yukawa plus linear potential v(r) = 10exp(— 5r)/r — 5exp(— r)/r + r. 

4 Numerical illustrations 

4.1 Yukawa plus linear confinement 

We illustrate the efficiency of our method with the example of non-relativistic kinetic 
energy with Yukawa plus linear confinement potential, 

v(r) — air + v\ exp(— \i\r)jr + V2 exp(—^2 r )/ r - (31) 

We choose the parameters as m = 1, a\ = 1, v\ = 10, pii = 5, v% = —5, and jj,2 = 1. 
This potential is shown in Fig. [l] Similar types of potentials are used to model the 
quark-quark interaction. We found that the parameters for the basis of the separable 
expansion 6 = 4 and b' = 3 were optimal. 

Table 2 The first three eigenstates in a Yukawa plus linear confinement in terms of basis 
states. 





Non-Relativistic 






Semi- 


Relativistic 




N 


1st 


2nd 


3rd 


N 


1st 


2nd 


3rd 


6 


-0.4380704 


2.396969 


3.820619 


18 


-1.596369 


1.007886 


2.324412 


9 


-0.4380688 


2.399815 


3.839923 


21 


-1.596349 


1.007905 


2.324431 


12 


-0.4380692 


2.399815 


3.840087 
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2.399815 
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2.399815 


3.840087 
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2.399815 
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-1.596355 
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2.324394 
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What is covered by 2 x 2 matrices is also covered by 3 x 3 matrices, therefore we 
use 3x3 blocks throughout. For a direct comparison of relativistic effects, we consider 
the non-relativistic and semi-relativistic Hamiltonian with the same potential. The 
convergence of the first three eigenstates of the non-relativistic and semi-relativistic 
Hamiltonian with increasing TV are shown in Table [2] To achieve a similar accuracy we 
need more basis states in the semi-relativistic case, but in both cases excellent accuracy 
has been observed by using a relatively small basis size. 



4.2 Yukawa plus quadratic confinement 

Another example we consider is Yukawa plus quadratic confinement potential (Fig. [2| 

v(r) = a 2 r 2 + v\ exp(— mr)/r + v 2 exp(—fi2r)/r , (32) 

with a non-relativistic and semi-relativistic kinetic energy operator. We use the param- 
eters m = 1, a<i = 0.5, i;i = 3, p,\ = 5, vi = —2, and \i% = \ , 




-3 J 

Fig. 2 Yukawa plus quadratic potential v(r) = 3exp(— 5r)/r — 2exp(— r)/r + 0.5r 2 . 

The convergence of the first three eigenstates of the non-relativistic and semi- 
relativistic Hamiltonian with increasing N are shown in Table [3] We can draw the 
same conclusion as before; excellent accuracy can be achieved using a a relatively 
small basis, but the convergence is not as rapid in the semi-relativistic case. 
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Table 3 The first three eigenstatcs of a Yukawa plus quadratic confining potential in terms 
of basis states. 



A? 



1st 



Non-Rclativistic 
2nd 3rd 



TV 



1st 



Semi-Rclativistic 
2nd 



3rd 



0.5502992 
0.5503106 
0.5503106 
0.5503106 
0.5503106 
0.5503106 



2.929190 
2.929148 
2.929150 
2.929150 
2.929150 
2.929150 



5.040812 
5.056162 
5.056224 
5.056225 
5.056225 
5.056225 



18 
21 
24 
27 
30 
33 



-0.1080278 
-0.1080290 
-0.1080297 
-0.1080301 
-0.1080304 
-0.1080306 



1.603202 
1.603201 
1.603200 
1.603200 
1.603200 
1.603199 



2.854743 
2.854741 
2.854741 
2.854740 
2.854740 
2.854740 



5 Summary and conclusions 



We developed a method for solving two-body non-relativistic and semi-relativistic prob- 
lems with confining potentials. We considered both linear and quadratic confinements. 
This method is based on a separable expansion of the short-range interaction in terms 
of CS basis and involves the evaluation of the corresponding Green's operator on that 
basis. The CS basis is particularly advantageous as it has simple analytic forms in 
configuration and momentum spaces allowing an easy evaluation of the potential and 
kinetic energy matrix elements. The non-relativistic kinetic energy plus confinement 
potential is represented by a infinite symmetric band matrix. The evaluation of the 
Green's operator amounts to inverting the infinite band matrix by means of matrix 
continued fractions. We found numerically that the semi-relativistic kinetic energy op- 
erator is asymptotically pentadiagonal in the CS basis. Therefore the Green's matrix 
of the semi-relativistic problem can also be evaluated using matrix continued fractions. 
As numerical examples we took the Yukawa plus linear and quadratic confinement 
models. We achieved very accurate results using a small number of basis functions. 
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